Single-cell and spatiotemporal transcriptomic analyses reveal the effects of microorganisms on immunity and metabolism in the mouse liver.

The gut-liver axis is a complex bidirectional communication pathway between the intestine and the liver in which microorganisms and their metabolites flow from the intestine through the portal vein to the liver and influence liver function. In a sterile environment, the phenotype or function of the liver is altered, but few studies have investigated the specific cellular and molecular effects of microorganisms on the liver. To this end, we constructed single-cell and spatial transcriptomic (ST) profiles of germ-free (GF) and specific-pathogen-free (SPF) mouse livers. Single-cell RNA sequencing (scRNA-seq) and single-nucleus RNA sequencing (snRNA-seq) revealed that the ratio of most immune cells was altered in the liver of GF mice; in particular, natural killer T (NKT) cells, IgA plasma cells (IgAs) and Kupffer cells (KCs) were significantly reduced in GF mice. Spatial enhanced resolution omics sequencing (Stereo-seq) confirmed that microorganisms mediated the accumulation of Kupffer cells in the periportal zone. Unexpectedly, IgA plasma cells were more numerous and concentrated in the periportal vein in liver sections from SPF mice but less numerous and scattered in GF mice. ST technology also enables the precise zonation of liver lobules into eight layers and three patterns based on the gene expression level in each layer, allowing us to further investigate the effects of microbes on gene zonation patterns and functions. Furthermore, untargeted metabolism experiments of the liver revealed that the propionic acid levels were significantly lower in GF mice, and this reduction may be related to the control of genes involved in bile acid and fatty acid metabolism. In conclusion, the combination of sc/snRNA-seq, Stereo-seq, and untargeted metabolomics revealed immune system defects as well as altered bile acid and lipid metabolic processes at the single-cell and spatial levels in the livers of GF mice. This study will be of great value for understanding host-microbiota interactions.

Microbiota-mediated shaping of mouse spleen structure and immune function characterized by scRNA-seq and Stereo-seq.

Gut microbes exhibit complex interactions with their hosts and shape an organism's immune system throughout its lifespan. As the largest secondary lymphoid organ, the spleen has a wide range of immunological functions. To explore the role of microbiota in regulating and shaping the spleen, we employ scRNA-seq and Stereo-seq technologies based on germ-free (GF) mice to detect differences in tissue size, anatomical structure, cell types, functions, and spatial molecular characteristics. We identify 18 cell types, 9 subtypes of T cells, and 7 subtypes of B cells. Gene differential expression analysis reveals that the absence of microorganisms results in alterations in erythropoiesis within the red pulp region and congenital immune deficiency in the white pulp region. Stereo-seq results demonstrate a clear hierarchy of immune cells in the spleen, including marginal zone (MZ) macrophages, MZ B cells, follicular B cells and T cells, distributed in a well-defined pattern from outside to inside. However, this hierarchical structure is disturbed in GF mice. Ccr7 and Cxcl13 chemokines are specifically expressed in the spatial locations of T cells and B cells, respectively. We speculate that the microbiota may mediate the structural composition or partitioning of spleen immune cells by modulating the expression levels of chemokines.

The first high-quality chromosome-level genome of the Sipuncula Sipunculus nudus using HiFi and Hi-C data.

Sipuncula is a class of exocoelomic unsegmented animals whose evolutionary relationships are unresolved. The peanut worm Sipunculus nudus is a globally distributed, economically important species belonging to the class Sipuncula. Herein, we present the first high-quality chromosome-level assembly of S. nudus based on HiFi reads and high-resolution chromosome conformation capture (Hi-C) data. The assembled genome was 1,427 Mb, with a contig N50 length of 29.46 Mb and scaffold N50 length of 80.87 Mb. Approximately 97.91% of the genome sequence was anchored to 17 chromosomes. A BUSCO assessment showed that 97.7% of the expectedly conserved genes were present in the genome assembly. The genome was composed of 47.91% repetitive sequences, and 28,749 protein-coding genes were predicted. A phylogenetic tree demonstrated that Sipuncula belongs to Annelida and diverged from the common ancestor of Polychaeta. The high-quality chromosome-level genome of S. nudus will serve as a valuable reference for studies of the genetic diversity and evolution of Lophotrochozoa.

Low-frequency and rare coding variants of NUS1 contribute to susceptibility and phenotype of Parkinson's disease.

NUS1 has been recently identified as a candidate gene for Parkinson's disease (PD). Few studies have examined the association of NUS1 variants with PD susceptibility and phenotypes. In the first cohort, whole-exome sequencing was performed to identify variants in NUS1 exon-coding and exon-intron regions in 1542 cases and 1625 controls. 13 variants were totally detected, of which 10 rare variants and 3 low-frequency variants. Burden analysis showed that rare NUS1 variants significantly enriched in PD (p=0.016). We also performed a meta-analysis based on previous and our studies to correlate NUS1 mutations with PD susceptibility. Integrating our previous cohort (3210 cases and 2807 controls) and the first cohort identified the significant association of rs539668656 with PD risk (odds ratio (OR) = 2.82, p = 0.016). The genotype-phenotype association analysis showed that patients carrying rare variants, or rs539668656 were significantly associated with earlier onset age, depression, emotional impairment and severe disease condition. Our results support the role of NUS1 rare variants and rs539668656 towards PD susceptibility and phenotype.

Nutrigenomics reveals potential genetic underpinning of diverse taste preference of Chinese men.

BACKGROUND: Taste preference varies geographically in China. However, studies on Chinese people's taste preference in different regions of China are limited, and are lack of research on the mechanism of differences in taste preference, especially in genetics. OBJECTIVE: This study aims to investigate the characteristics of taste preference of Chinese men, and estimate whether diverse taste preference in Chinese have genetic underpinning. METHODS: We conducted a questionnaire survey on taste preferences on 1076 males from 10 regions of China, and collected another 1427 males from the same regions which genotyped by microarray. We compared the correlation between different taste preference, and evaluated the correlation between the mutation frequency of inhouse database and different taste preference. The putative taste-preference-related genes were further utilized to estimate the candidate relationship on gene and gene network in different taste preference. RESULTS: There was a correlation between different taste preferences in Chinese men. We found 31 SNPs associated with 6 kind of taste preferences. These SNPs located within or nearby 36 genes, and the tastes associated with 4 of these genes (TRPV1, AGT, ASIC2 and GLP1R) are consistent with the previous studies. Moreover, in different tastes which were suggested to be associated with each other, some putative related genes were the same or in the same gene network, such as pathways related with blood pressure, response to stimulus and nervous system. CONCLUSIONS: This study indicates that the diverse taste preference of Chinese men may have genetic underpinning.

Genome-Wide Association Study of Smoking Behavior Traits in a Chinese Han Population.

Tobacco use is one of the leading causes of preventable disease worldwide. Genetic studies have elucidated numerous smoking-associated risk loci in American and European populations. However, genetic determinants for cigarette smoking in Chinese populations are under investigated. In this study, a whole-genome sequencing (WGS)-based genome-wide association study (GWAS) was performed in a Chinese Han population comprising 620 smokers and 564 nonsmokers. Thirteen single-nucleotide polymorphisms (SNPs) of the raftlin lipid linker 1 (RFTN1) gene achieved genome-wide significance levels (P < 5 x 10-8) for smoking initiation. The rs139753473 from RFTN1 and six other suggestively significant loci from CUB and sushi multiple domains 1 (CSMD1) gene were also associated with cigarettes per day (CPD) in an independent Chinese sample consisting of 1,329 subjects (805 smokers and 524 nonsmokers). When treating males separately, associations between smoking initiation and PCAT5/ANKRD30A, two genes involved in cancer development, were identified and replicated. Within RFTN1, two haplotypes (i.e., C-A-C-G and A-G-T-C) formed by rs796812630-rs796584733-rs796349027-rs879511366 and three haplotypes (i.e., T-T-C-C-C, T-T-A-T-T, and C-A-A-T-T) formed by rs879401109-rs879453873-rs75180423-rs541378415-rs796757175 were strongly associated with smoking initiation. In addition, we also revealed two haplotypes (i.e., C-A-G-G and T-C-T-T derived from rs4875371-rs4875372-rs17070935-rs11991366) in the CSMD1 gene showing a significant association with smoking initiation. Further bioinformatics functional assessment suggested that RFTN1 may participate in smoking behavior through modulating immune responses or interactions with the glucocorticoid receptor alpha and the androgen receptor. Together, our results may help understand the mechanisms underlying smoking behavior in the Chinese Han population.

GCH1 variants contribute to the risk and earlier age-at-onset of Parkinson's disease: a two-cohort case-control study.

BACKGROUND: Common and rare variants of guanosine triphosphate cyclohydrolase 1 (GCH1) gene may play important roles in Parkinson's disease (PD). However, there is a lack of comprehensive analysis of GCH1 genotypes, especially in non-coding regions. The aim of this study was to explore the genetic characteristics of GCH1, including rare and common variants in coding and non-coding regions, in a large population of PD patients in Chinese mainland, as well as the phenotypic characteristics of GCH1 variant carriers. METHODS: In the first cohort of this case-control study, we performed whole-exome sequencing in 1555 patients with early-onset or familial PD and 2234 healthy controls; then in the second cohort, whole-genome sequencing was performed in sporadic late-onset PD samples (1962 patients), as well as 1279 controls. Variants at target GCH1 regions were extracted, and then genetic and detailed phenotypic data were analyzed using regression models and the sequence kernel association test. We also performed a meta-analysis to correlate deleterious GCH1 variants with age at onset (AAO) in PD patients. RESULTS: For coding variants, we identified a significant burden of GCH1 deleterious variants in early-onset or familial PD cases compared to controls (1.2% vs 0.1%, P < 0.0001). In the analysis of possible regulatory variants in GCH1 non-coding regions, rs12323905 (P = 0.001, odds ratio = 1.19, 95%CI 1.07-1.32) was significantly associated with PD, and variant sets in untranslated regions and intron regions, GCH1 brain-specific expression quantitative trait loci, and two possible promoter/enhancer (GH14J054857 and GH14J054880) were suggestively associated with PD. Genotype-phenotype correlation analysis revealed that the carriers of GCH1 deleterious variants manifested younger AAO (P < 0.0001), and had milder motor symptoms, milder fatigue symptoms and more autonomic nervous dysfunctions. Meta-analysis of six studies demonstrated 6.4-year earlier onset in GCH1 deleterious variant carriers (P = 0.0009). CONCLUSIONS: The results highlight the importance of deleterious variants and non-coding variants of GCH1 in PD in Chinese mainland and suggest that GCH1 mutation can influence the PD phenotype, which may help design experimental studies to elucidate the mechanisms of GCH1 in the pathogenesis of PD.

High-throughput sequencing reveals the diversity of TCR ß chain CDR3 repertoire in patients with severe acne.

Acne is a common chronic inflammatory skin disease, and the inflammation immune response runs through all stages of acne lesions. In this study, we use a combination of multiplex-PCR and high-throughput sequencing technologies to analyze T cell receptor ß chain CDR3 (complementarity-determining region 3) in peripheral blood isolated from severe acne patients. Once compared with healthy controls, we propose to identify acne-relevant CDR3 peptides. Our results reveal that the diversity of T cell receptor ß chain (TRB) CDR3 sequences in the peripheral blood of the severe acne vulgaris (SA) group differed from that of the control group. In addition, we find 10 TRB CDR3 sequences, amino acid sequences and V-J combinations with significantly different expressions between the SA group and the non-acne (NA) group (P < 0.0001). These findings may contribute to a better understanding of the role of immunity in the pathogenesis of acne and may serve as biomarkers for evaluating risk or prognosis of severe acne disease in future.

Genomic landscapes of Chinese sporadic autism spectrum disorders revealed by whole-genome sequencing.

Autism spectrum disorder (ASD) is a neurodevelopmental disorder with considerable clinical and genetic heterogeneity. In this study, we identified all classes of genomic variants from whole-genome sequencing (WGS) dataset of 32 Chinese trios with ASD, including de novo mutations, inherited variants, copy number variants (CNVs) and genomic structural variants. A higher mutation rate (Poisson test, P < 2.2 × 10-16) in exonic (1.37 × 10-8) and 3'-UTR regions (1.42 × 10-8) was revealed in comparison with that of whole genome (1.05 × 10-8). Using an integrated model, we identified 87 potentially risk genes (P < 0.01) from 4832 genes harboring various rare deleterious variants, including CHD8 and NRXN2, implying that the disorders may be in favor to multiple-hit. In particular, frequent rare inherited mutations of several microcephaly-associated genes (ASPM, WDR62, and ZNF335) were found in ASD. In chromosomal structure analyses, we found four de novo CNVs and one de novo chromosomal rearrangement event, including a de novo duplication of UBE3A-containing region at 15q11.2-q13.1, which causes Angelman syndrome and microcephaly, and a disrupted TNR due to de novo chromosomal translocation t(1; 5)(q25.1; q33.2). Taken together, our results suggest that abnormalities of centrosomal function and chromatin remodeling of the microcephaly-associated genes may be implicated in pathogenesis of ASD. Adoption of WGS as a new yet efficient technique to illustrate the full genetic spectrum in complex disorders, such as ASD, could provide novel insights into pathogenesis, diagnosis and treatment.

Coding mutations in NUS1 contribute to Parkinson's disease.

Whole-exome sequencing has been successful in identifying genetic factors contributing to familial or sporadic Parkinson's disease (PD). However, this approach has not been applied to explore the impact of de novo mutations on PD pathogenesis. Here, we sequenced the exomes of 39 early onset patients, their parents, and 20 unaffected siblings to investigate the effects of de novo mutations on PD. We identified 12 genes with de novo mutations (MAD1L1, NUP98, PPP2CB, PKMYT1, TRIM24, CEP131, CTTNBP2, NUS1, SMPD3, MGRN1, IFI35, and RUSC2), which could be functionally relevant to PD pathogenesis. Further analyses of two independent case-control cohorts (1,852 patients and 1,565 controls in one cohort and 3,237 patients and 2,858 controls in the other) revealed that NUS1 harbors significantly more rare nonsynonymous variants (P = 1.01E-5, odds ratio = 11.3) in PD patients than in controls. Functional studies in Drosophila demonstrated that the loss of NUS1 could reduce the climbing ability, dopamine level, and number of dopaminergic neurons in 30-day-old flies and could induce apoptosis in fly brain. Together, our data suggest that de novo mutations could contribute to early onset PD pathogenesis and identify NUS1 as a candidate gene for PD.

Identification of MicroRNA Targets of Capsicum spp. Using MiRTrans-a Trans-Omics Approach.

The microRNA (miRNA) can regulate the transcripts that are involved in eukaryotic cell proliferation, differentiation, and metabolism. Especially for plants, our understanding of miRNA targets, is still limited. Early attempts of prediction on sequence alignments have been plagued by enormous false positives. It is helpful to improve target prediction specificity by incorporating the other data sources such as the dependency between miRNA and transcript expression or even cleaved transcripts by miRNA regulations, which are referred to as trans-omics data. In this paper, we developed MiRTrans (Prediction of MiRNA targets by Trans-omics data) to explore miRNA targets by incorporating miRNA sequencing, transcriptome sequencing, and degradome sequencing. MiRTrans consisted of three major steps. First, the target transcripts of miRNAs were predicted by scrutinizing their sequence characteristics and collected as an initial potential targets pool. Second, false positive targets were eliminated if the expression of miRNA and its targets were weakly correlated by lasso regression. Third, degradome sequencing was utilized to capture the miRNA targets by examining the cleaved transcripts that regulated by miRNAs. Finally, the predicted targets from the second and third step were combined by Fisher's combination test. MiRTrans was applied to identify the miRNA targets for Capsicum spp. (i.e., pepper). It can generate more functional miRNA targets than sequence-based predictions by evaluating functional enrichment. MiRTrans identified 58 miRNA-transcript pairs with high confidence from 18 miRNA families conserved in eudicots. Most of these targets were transcription factors; this lent support to the role of miRNA as key regulator in pepper. To our best knowledge, this work is the first attempt to investigate the miRNA targets of pepper, as well as their regulatory networks. Surprisingly, only a small proportion of miRNA-transcript pairs were shared between degradome sequencing and expression dependency predictions, suggesting that miRNA targets predicted by a single technology alone may be prone to report false negatives.

PEAR1 gene polymorphism in a Chinese pedigree with pulmonary thromboembolism.

To explore the correlation between platelet endothelial aggregation receptor-1 (PEAR1) genetic polymorphism and pulmonary thromboembolism (PTE).Variant loci of the PEAR1 gene were screened in a PTE pedigree, followed by verification using Sanger sequencing. These polymorphic loci were validated in 101 PTE patients and 132 matched normal patients using MassARRAY single nucleotide polymorphism (SNP) genotyping methods. The frequency differences between the allele and genotypes were compared using the Hardy-Weinberg equilibrium test and Chi-square test. The correlation between the PEAR1 gene SNP and PTE was analyzed by comparing the between-group variance differences using the χ test.Three SNPs were identified in the PTE pedigree. There was a heterozygous transition of T>C in rs1952294, and a transition of C>T in rs778026543 in 2 members in the pedigree; however, the rs778026543 was not identified in the 101 PTE patients and 132 healthy controls. The genotype and allele frequencies of rs822442 did not differ significantly between PTE patients and healthy controls (P > 0.05). The variance difference at rs778026543 between pedigree members and healthy controls was significant (P < 0.001), supporting its potential heredity.The PEAR1 polymorphism, rs778026543, but not rs1952294 and rs822442, may be a susceptibility SNP for PTE.

Deep sequencing of the MHC region in the Chinese population contributes to studies of complex disease.

The human major histocompatibility complex (MHC) region has been shown to be associated with numerous diseases. However, it remains a challenge to pinpoint the causal variants for these associations because of the extreme complexity of the region. We thus sequenced the entire 5-Mb MHC region in 20,635 individuals of Han Chinese ancestry (10,689 controls and 9,946 patients with psoriasis) and constructed a Han-MHC database that includes both variants and HLA gene typing results of high accuracy. We further identified multiple independent new susceptibility loci in HLA-C, HLA-B, HLA-DPB1 and BTNL2 and an intergenic variant, rs118179173, associated with psoriasis and confirmed the well-established risk allele HLA-C*06:02. We anticipate that our Han-MHC reference panel built by deep sequencing of a large number of samples will serve as a useful tool for investigating the role of the MHC region in a variety of diseases and thus advance understanding of the pathogenesis of these disorders.

Whole-exome sequencing implicates UBE3D in age-related macular degeneration in East Asian populations.

Age-related macular degeneration (AMD) is a leading cause of irreversible central blindness among the elderly worldwide. We use exome sequencing to analyse nonsynonymous single-nucleotide variants (SNVs) across the whole genome of 216 neovascular AMD cases and 1,553 controls. As a follow-up validation, we evaluate 3,772 neovascular AMD cases and 6,942 controls from five independent cohorts in the East Asian population. Here we show strong evidence of an association at a novel, missense SNV, rs7739323, which is located in the ubiquitin protein ligase E3D (UBE3D) gene (Pmeta=1.46 × 10(-9), odds ratio (OR)=0.74, 95% confidence interval (CI): 0.63-0.88). Furthermore, ablation of the UBE3D protein lead to an abnormal amount of pigment granules deposited in retinal pigment epithelium microvilli area and an abnormal response on electroretinography (ERG) in UBE3D(+/-) heterozygous mice. Our findings indicate that the ubiquitin-proteasome system may play a role in the pathogenesis of neovascular AMD.

Genetic diagnosis of two dopa-responsive dystonia families by exome sequencing.

Dopa-responsive dystonia, a rare disorder typically presenting in early childhood with lower limb dystonia and gait abnormality, responds well to levodopa. However, it is often misdiagnosed with the wide spectrum of phenotypes. By exome sequencing, we make a rapid genetic diagnosis for two atypical dopa-responsive dystonia pedigrees. One pedigree, presented with prominent parkinsonism, was misdiagnosed as Parkinson's disease until a known mutation in GCH1 (GTP cyclohydrolase 1) gene (NM_000161.2: c.631_632delAT, p.Met211ValfsX38) was found. The other pedigree was detected with a new compound heterozygous mutation in TH (tyrosine hydroxylase) gene [(NM_000360.3: c.911C>T, p.Ala304Val) and (NM_000360.3: c.1358G>A, p.Arg453His)], whose proband, a pregnant woman, required a rapid and less-biased genetic diagnosis. In conclusion, we demonstrated that exome sequencing could provide a precise and rapid genetic testing in the diagnosis of Mendelian diseases, especially for diseases with wide phenotypes.

Sequencing-based approach identified three new susceptibility loci for psoriasis.

In a previous large-scale exome sequencing analysis for psoriasis, we discovered seven common and low-frequency missense variants within six genes with genome-wide significance. Here we describe an in-depth analysis of noncoding variants based on sequencing data (10,727 cases and 10,582 controls) with replication in an independent cohort of Han Chinese individuals consisting of 4,480 cases and 6,521 controls to identify additional psoriasis susceptibility loci. We confirmed four known psoriasis susceptibility loci (IL12B, IFIH1, ERAP1 and RNF114; 2.30 × 10(-20)≤P≤2.41 × 10(-7)) and identified three new susceptibility loci: 4q24 (NFKB1) at rs1020760 (P=2.19 × 10(-8)), 12p13.3 (CD27-LAG3) at rs758739 (P=4.08 × 10(-8)) and 17q12 (IKZF3) at rs10852936 (P=1.96 × 10(-8)). Two suggestive loci, 3p21.31 and 17q25, are also identified with P<1.00 × 10(-6). The results of this study increase the number of confirmed psoriasis risk loci and provide novel insight into the pathogenesis of psoriasis.

A large-scale screen for coding variants predisposing to psoriasis.

To explore the contribution of functional coding variants to psoriasis, we analyzed nonsynonymous single-nucleotide variants (SNVs) across the genome by exome sequencing in 781 psoriasis cases and 676 controls and through follow-up validation in 1,326 candidate genes by targeted sequencing in 9,946 psoriasis cases and 9,906 controls from the Chinese population. We discovered two independent missense SNVs in IL23R and GJB2 of low frequency and five common missense SNVs in LCE3D, ERAP1, CARD14 and ZNF816A associated with psoriasis at genome-wide significance. Rare missense SNVs in FUT2 and TARBP1 were also observed with suggestive evidence of association. Single-variant and gene-based association analyses of nonsynonymous SNVs did not identify newly associated genes for psoriasis in the regions subjected to targeted resequencing. This suggests that coding variants in the 1,326 targeted genes contribute only a limited fraction of the overall genetic risk for psoriasis.

Detection of clinically relevant genetic variants in autism spectrum disorder by whole-genome sequencing.

Autism Spectrum Disorder (ASD) demonstrates high heritability and familial clustering, yet the genetic causes remain only partially understood as a result of extensive clinical and genomic heterogeneity. Whole-genome sequencing (WGS) shows promise as a tool for identifying ASD risk genes as well as unreported mutations in known loci, but an assessment of its full utility in an ASD group has not been performed. We used WGS to examine 32 families with ASD to detect de novo or rare inherited genetic variants predicted to be deleterious (loss-of-function and damaging missense mutations). Among ASD probands, we identified deleterious de novo mutations in six of 32 (19%) families and X-linked or autosomal inherited alterations in ten of 32 (31%) families (some had combinations of mutations). The proportion of families identified with such putative mutations was larger than has been previously reported; this yield was in part due to the comprehensive and uniform coverage afforded by WGS. Deleterious variants were found in four unrecognized, nine known, and eight candidate ASD risk genes. Examples include CAPRIN1 and AFF2 (both linked to FMR1, which is involved in fragile X syndrome), VIP (involved in social-cognitive deficits), and other genes such as SCN2A and KCNQ2 (linked to epilepsy), NRXN1, and CHD7, which causes ASD-associated CHARGE syndrome. Taken together, these results suggest that WGS and thorough bioinformatic analyses for de novo and rare inherited mutations will improve the detection of genetic variants likely to be associated with ASD or its accompanying clinical symptoms.